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Abstract. Based on XMM-Newton EPIC data of four pointings towards the Small Magellanic Cloud (SMC), 
results on timing and spectral analyses of 16 known high mass X-ray binaries (HMXBs) and HMXB candidates 
in the SMC are presented. We confirm the pulse periods of four sources which were known to show pulsations. 
In addition, two new X-ray pulsars are discovered: XMMU J005605. 2- 722200 with P pu i 8e = 140.1 ± 0.3 s and 
RX J0057. 8-7207 with P pu i sc = 152.34 ± 0.05 s. Due to the low Galactic foreground absorption, X-ray binary 
systems in the Magellanic Clouds are well suited for studies of the soft component in their X-ray spectrum. Spectral 
analysis reveals soft emission besides a power law component in the spectra of three sources. The existence of 
emission lines in at least one of them corroborates the thermal nature of this emission with temperatures of 0.2 - 
0.3 keV and heavy element abundances lower than solar. For the HMXB SMCX-2 which was in a low luminosity 
state, we determine a flux upper limit of 1.5 x 10 -14 erg cm -2 s _1 (0.3 - 10.0 keV). Furthermore, two new sources 
(XMMUJ005735. 7-721932 and XMMU J010030. 2-722035) with hard spectrum and emission line objects as likely 
optical counterparts are proposed as new X-ray binary candidates. 
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1. Introduction 

After the discovery of X-ra y emission f r om the Magellanic 
Clouds (MCs) in 1970 (Price et all Il97l[) . surveying 
observations of each MC were performed by differ- 
ent X-ray observatories. As for the Small Magellanic 
Cloud (SMC), source catal ogues were created from 



observations wit h Eins t ein llSeward fc Mitchell Il98ll 



(Kahabka et all Il999: 


Haberl et al 


. 1200(1 ISasaki et all 


120001). and ASCA dYokosrawa et all 


2000). 



The analysis of these X-ray sources has shown, that a 
large number of X-ray bright objects belongs to the class 
of X-ray binaries (XRBs) in which a neutron star or a 
black hole forms a binary system with a companion star. 
In these systems, mass is accreted from the donor star 
onto the compact object. X-ray binaries can be divided 
into low mass X-ray binaries and high mass (or massive) 
X-ray binaries (HMXBs), depending on the mass of the 
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companion star. Therefore, the identification of optical 
counterparts of the X-ray sources is crucial for the un- 
derstanding of the nature of these sources. Furthermore, 
HMXBs form two subgroups with either an OB supergiant 
or a Be star as donor. A detail e d catalogue of HMXBs was 
compiled by iLiu et all l)2000|) . iNeeueruela k Cod l)2002l) 
performed high resolution spectroscopy of optical coun- 
terparts of HMXBs in the Large Magellanic Cloud (LMC) 
and studied the population of HMXBs. In the Milky Way 
or in the LMC, the fraction of Be/X-ray binary systems 
(Be/XRB) is 60 - 70% of all HMXBs, whereas more than 
90% of the HMXBs in the SM C turned out to be Be sys- 
tems IH aberl fc Sasaki l2000l and references therein). 

Since pulsed X-ray emission can be observed from 
neutron star HMXBs, these sources are also called X- 
ray binary pulsars. Based on ASCA, RXTE, ROSAT 
and Beppo SAX observations, more than 20 X-ray bi- 
nary pulsars have been discovered in t he SM C so far 
l|Haberl fc Sasaki! l200fj Fibkogawa et"afl 1200(1 and refer- 
ences therein) . Moreover, in one of t he first observations of 
XMM-Newton l|Jansen et all 120011) , pulsed emission from 
anot her HMXB was foun d, which was identified with a Be 
star ((Sasaki et alll200lh . 

In order to improve our understanding of the X-ray 
source population in the SMC, we proposed and analy- 
sed pointed observations of the SMC by XMM-Newton 
and performed spectral and temporal studies of detected 
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Fig. 1. DSS image of the SMC and the position of the 
EPIC field of view of the XMM-Newton observations listed 
in Tabled The discontinuity seen in the DSS image is an 
artifact. 



sources. In this paper, we focus on the class of HMXBs 
and present the results on each HMXB and candidate in 
the observed fields. 



2. Data 

For AO-1 of XMM-Newton, we proposed observations 
of eight fields in the SMC in order to study the X-ray 
binary population (PI: W.P.). Two observations of this 
proposal were performed. During the first observation 
(ID 00842008), the telescope pointed towards the HMXB 
SMC X-2 in the south of the main galaxy. In order to per- 
form source detection and analysis in the whole field of 
view, we used data from the European P hoton Imaging 
Cameras EPIC P N llStriider et all 12001ft . EPIC MOS1, 
and EPIC MOS2 l|Turner et all l2Q0lh . The observation 
was performed with all the EPIC cameras in full frame 
mode. For EPIC PN the thin filter was used, whereas for 
the EPIC MOS cameras medium filters were chosen. The 
next observation, ID 00842001, covered a region in the 
north of the SMC. The CCD read out modes and the fil- 
ters of the EPIC cameras were the same as in the first 
observation. 

Moreover, we searched the XMM-Newton Science Data 
Archive for public data of the SMC, suitable for our 
purposes. We found two data sets (ID 01100002 and 
01357206) of fields in the north of the SMC, slightly over- 
lapping with each other as well as with the pointing ID 
00842001. The details of the observations are summarised 
in Tabled 



2.1. Source detection 

All the data were processed with the XMM-Newton 
Science Analysis System (SAS) version 5.3.3. For source 



detection, the events were separated into four energy 
bands: B 1 = 0.3 - 1.0 keV, B 2 = 1.0 - 2.0 keV, B 3 = 2.0 - 
4.5 keV, and B 4 = 4.5 - 10.0 keV. In all these bands, im- 
ages were created and source detection was performed us- 
ing the sliding window and maximum likelihood methods 
of the SAS. Detections with likelihood of existence (ML) 
higher than 10.0 were accepted as real sources. This corre- 
sponds to the probability P = 1 - cxp(-ML) = 0.999955 
for the existence of the source. Hardness ratios were com- 
puted using the source counts in different bands: 



HRi = 



B H 



B, 



B i+ i + Bi 



(1) 



for i = 1,2,3. The values of HRI, HR2, and HR3 are 
shown in two diagrams in Fig. [2] In most cases, X-ray bi- 
naries in the SMC or Active Galactic Nuclei (AGNs) be- 
hind the SMC have absorbed spectra and therefore show 
positive values for HRI. As can be seen in the second di- 
agram, HR2 and HR3 which compare the events in the 
energy bands above 1.0 keV, cluster around zero. A large 
fraction of the sources (90%) is located in a region with 
-0.4 < HR2 < 0.3 and -0.6 < HR3 < 0.5. 

The detected sources w ere cross-correlated with cata- 
logues of Einstein sources tWang fc Wul.ll992i) as well as 
RQSAT sour c es det ected bv PSP C faahabka et allll99flt 
lHaberl et all l2000h and by HRI <|Sasaki et all l2000h in- 
struments. The positions of the X-ray sources were plotted 
on Digitized Sky Survey DSS2 (red) images of this field 
in order to find probable optical counterparts. The opti- 
cal sources were also verified by cross-correlating the X-ray 
source list with the USNO-A2.0 ca talogue produced by the 
United States Naval Observatory l|Monetl ll 99*1 1 1998(1 . In 
addition, we compared the source positions to the entries 
in the Optical Gravitational Lensing Experiment OGLE- 
II project list of varia ble sources in the Magellanic Clouds 
l(Zebrun et all l200l(l . For five out of 15 sources, corre- 
lations with an OGLE object were found. Finally, the 
source list was cross-correlated with the list of e mission 
line objects in the SMC (Meyssonnier & Azzopardi. Il993l 
[MA93] ) . The existence of an emission line star at the posi- 
tion of a hard X-ray source indicates that the source might 
be an X-ray binary system with a Be-star companion. 

The complete set of source lists will be presented in 
another paper. Here, we shall concentrate on the HMXBs 
and candidates in the four fields. The results on the eigh- 
teen sources are summarised in Table The table in- 
cludes the X-ray source coordinates, 1 a positional er- 
ror, count rate, likelihood for detection in the total band, 
flux, hardness ratios (see Eq. Q), pulse period with la- 
error fSect. l2~2l . and identifications. Correlations with 
RQSAT sources, OGLE objects, a nd emission line objects 
l|Mev ssonnier fc Azz opardil 1199^ can be found as well. 
The sources are sorted by RA and Dec (J2000.0), and the 
entry numbers are used in the following. The positional 
errors which are given in this table are statistical errors. 
The systematic error of t he X-ray position is about 3" 
- 4" IfBarcons et al ., 2002). In order to calculate the flux, 
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Table 1. XMM-Newton data used for the analysis. 



1 


2 


3 


4 


5 


6 


7 


8 


Rev. 


Obs. ID 


Pointing direction 


Inst. 


Mode 


Filter 


Start time (UT) 


End time (UT) 






RA Dec 


(EPIC) 














(J2000.0) 












157 


01100002 


00 59 46.6 -72 09 30 


PN 


Full 


Medium 


2000/10/17 16:16:36 


2000/10/17 20:41:09 








Ml 


Full 


Medium 


2000/10/17 15:10:44 


2000/10/17 20:39:43 








M2 


Full 


Medium 


2000/10/17 15:10:35 


2000/10/17 20:39:42 


247 


01357206 


01 03 29.0 -72 02 33 


PN 


Full 


Thinl 


2001/04/15 01:20:28 


2001/04/15 05:50:27 








Ml 


PW3 


Thinl 


2001/04/14 20:47:25 


2001/04/15 05:55:45 








M2 


PW3 


Thinl 


2001/04/14 20:47:25 


2001/04/15 05:55:45 


340 


00842008 


00 54 54.3 -73 40 12 


PN 


Full 


Thinl 


2001/10/17 10:46:50 


2001/10/17 15:55:11 








Ml 


Full 


Medium 


2001/10/17 10:07:40 


2001/10/17 15:59:35 








M2 


Full 


Medium 


2001/10/17 10:07:40 


2001/10/17 15:59:36 


422 


00842001 


00 56 24.4 -72 21 33 


PN 


Full 


Thinl 


2002/03/30 14:21:45 


2002/03/30 19:39:38 








Ml 


Full 


Medium 


2002/03/30 13:48:28 


2002/03/30 19:44:36 








M2 


Full 


Medium 


2002/03/30 13:48:29 


2002/03/30 19:44:54 



Notes to column No 5: PW3: Partial Window 3. 
Notes to column No 6: Ml: MOS1, M2: MOS2. 



model parameters resulting from the spectral analysis (see 
Sect. 12 .21) were used for all sources. 

In the following, the s ource number in t he ROSAT 
HRI catalogue of the SMC l|Sasaki et all 12000^) is given as 
RHNNN, and the n umber in the ROSAT PSPC catalogue 
l|Haberl et alll2000f) as RPNNN. The entry number in the 
list of Haberl & Sasaki l|200(i [HS2000]) is also mentioned 
using the format [HS2000] NN (also see Table EJ 

2.2. Timing and spectral analysis 

As the very first step of data analysis, we checked the 
EPIC PN data for incorrect time information. It has been 
reported that in some cases, there are time jumps of 1 s in 
the EPIC PN data which were not corrected in the SAS 
processing. Since we didn't find any event which indicated 
such a time jump, we could proceed without any counter- 
measure. 

After selecting the events for each source, they were 
analysed using the XANADU software package dis- 
tributed by the High Energy Astrophysics Science Archive 
Research Center (HEASARC). It contains the packages 
XRONOS for timing analysis and XSPEC for spectral fit- 
ting. 

Based on EPIC PN data, period search was carried 
out with XRONOS after correcting the photon arrival 
times for solar system barycentre. If a peak was found in 
the power spectrum indicating pulsations, a more detailed 
epoch folding search was performed around the prelimi- 
nary value. Once we got the rough value for the pulse pe- 
riod, the x 2 distribution around this value was fitted with 
a Lorentz profile and the maximum of the Lorentz profile 
was determined together with the 1 tr error. Finally, folded 
light curves were created in three energy bands: Bi = 0.3 
- 1.0 keV, B 2 = 1.0 - 2.0 keV, B 3+4 = 2.0 - 10.0 keV. In 
addition, the ratio of the count rates in the harder band to 



the softer band (B 2 /Bi and B 3+4 /Bi +2 , with Bi +2 = 0.3 
- 2.0 keV) were computed to illustrate the changes in the 
hardness ratios with pulse phase. Note that these hardness 
ratios are different to the numbers defined in Eq. 

Except for sources which were too faint, spectra were 
extracted for each source. These spectra were modelled 
with a power law component together with the fixed 
Galactic foreground absorbing column density of iVH,Gai = 
5.74 x 10 20 cm -2 l)Dickev fc Lockmanl Il990|) and a free 
column density iVjj: 

S^E) = e-^^HG.i x e -°(E)N H xK x E -r^ ^ 

E being the energy in [keV], T the photon index, and 
K the normalisation. In some cases there was a devia- 
tion of the observed spectrum from a power law spectrum, 
suggesting the existence of an additional soft component. 
Since the spectra show features indicating emission lines, 
the soft component was modelled as thermal plasma emis- 
sion. This thermal emission presumably arises from cir- 
cumstellar matter, and the absorption must be negligibly 
low in comparison to the absorption of the hard X-ray 
emission from the neutron star. Moreover, if the column 
density was high, the soft emission would be absorbed and 
thus not detectable . Using the MEKAL model in XSPEC 
llMewe et all Il985l Il986t iKaastraL Il992t iLiedahl et all 
Il995h for the thermal component without additional ab- 
sorbing column density, the spectrum can be written as 

S 2 {E) = e-^^HGai x [ e -<m^ x R x E- T 

+Smekal(T, Abund.)). (3) 

•SmekalC^) i s the MEKAL model spectrum with a tem- 
perature corresponding to kT in [keV] and elemental 
abundances with respect to solar. We also performed a 
fit with a blackbody component instead of the MEKAL 
model. Although for fainter sources, no significant differ- 
ence was found in the fits, for bright sources like No 16 
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Table 2. Detected HMXBs and candidates in the fields of the SMC observed by XMM-Newton. 



1 

J- 


2 


3 


4 







6 




7 


« 





No 


Obs. ID 


RA 


Dec 


Pos. err. 


Count rate 




Det. ML 


Flux 








(J2000.0) 


["] 




Is" 1 ] 






[erg cm" 2 




01 


00842008 Ml 


00 51 56.05 


-73 41 51.4 


2.1 


6.33 x 10" 


3 ± 1.04 x 10" 


■3 


84.2 


~ 1 x 10" 


-13 


02 


00842008 PN 


00 54 33.4 


-73 41 04° 


- 


< 2.33 x 10~ 03t 




3.4 


< 1.5 x 10 


-14 


03 


00842001 PN 


00 54 56.02 


-72 26 48.6 


1.0 


3.30 x 10" 


2 ±3.44 x 10" 


-:i 


267.1 


5.3 x 10" 


14 


04 


00842001 PN 


00 56 05.24 


-72 22 00.9 


2.0 


7.46 x 10" 


3 ± 1.32 x 10" 


-:i 


76.0 


2.4 x 10" 


14 


05 


00842001 PN 


00 57 19.58 


-72 25 35.1 


0.5 


8.79 x 10" 


2 ±4.16 x 10" 


■3 


225.6 


1.4 x 10" 


13 


06 


00842001 PN 


00 57 35.71 


-72 19 32.6 


0.9 


3.07 x 10" 


2 ± 2.46 x 10" 


■3 


478.5 


9.5 x 10" 


14 


06a 


01100002 PN 


00 57 35.56 


-72 19 36.8 


2.9 


2.53 x 10" 


2 ±6.94 x 10" 


■3 


55.9 


1.9 x 10" 


14 


07 


01100002 PN 


00 57 50.22 


-72 02 37.0 


0.7 


1.32 x 10" 


1 ± 8.75 x 10" 


■3 


1147.0 


2.3 x 10" 


13 


08 


01100002 PN 


00 57 50.80 


-72 07 58.7 


0.5 


1. / y x iu 


■1 1 Q 11 1 n - 

zh o.ll X IU 


■3 


2798.1 


5.8 x 10" 


13 


09 


00842001 PN 


00 58 11.68 


-72 30 50.4 


1.7 


5.37 x 10" 


2 ±8.49 x 10" 


-3 


120.2 


~ 3 x 10" 


-13 


10 


01100002 Ml 


nn p;o 01 ni 
UU Dy zl.Ul 


70 OQ 1 C A 
— ( Z Zo lo.4 


U.O 


5.22 x 10" 


2 ±3.22 x 10" 


-:i 


1269.6 


z.o X IU 


13 


10a 


00842001 PN 


00 59 21.10 


-72 23 15.8 


0.9 


1 nc; \^ in - 
l.UO X iu 


"1-1-7 OQ \s 1 n~ 
it i .Zo X iu 


■3 


828.9 


1.2 x 10" 


13 


1 1 


ai i nnnni r>i\T 
UllUUUUz -rlN 


m nn 1 a 1 c 
Ul UU lo.lo 


70 f\A A ^ Q 
— ( Z U4 40. 


l.z 


/in \y 1 n ~ 
z.4z X IU 


zh Z.oO X IU 


■3 


A01.O 


a ^, 1 n - 
4.z X IU 


14 


12 


01100002 PN 


m nn oq 
ui uu ou.zo 


70 on qp; 1 

— iZ Z\J 00. 1 


3.2 


o 01 v- 1 n~ 
y.oi x iu 


3 _i_ q nn v in - 
zt o.uu X iu 


-3 


25.8 


~ 6 x 10" 




13 


01357206 PN 


ni ni no os 
ui ui uz.yo 


70 Ofi t^Q K 

— ( z ud oy.o 


2.9 


1.77 x 10" 


2 ±4.16 x 10" 


■3 


38.2 


2.6 x 10" 


14 


13a 


01100002 PN 


01 01 03.87 


-72 07 02.4 


3.8 


9.16 x 10" 


3 ± 2.68 x 10" 


-3 


23.0 


3.2 x 10" 




14 


01100002 PN 


01 01 20.82 


-72 11 21.1 


0.6 


1.40 x 10" 


1 ± 7.47 x 10" 




1689.6 


3.1 x 10" 




14a 


01357206 PN 


01 01 20.87 


-72 11 16.8 


1.0 


8.83 x 10" 


2 ±8.74 x 10" 


-3 


368.7 


8.5 x 10" 


14 


15 


01357206 PN 


01 01 37.56 


-72 04 18.7 


1.0 


5.01 x 10" 


2 ±4.26 x 10" 


■3 


436.3 


5.7 x 10" 


14 


15a 


01100002 PN 


01 01 37.77 


-72 04 22.1 


1.2 


3.55 x 10" 


2 ±4.36 x 10" 


■3 


222.7 


4.6 x 10" 


14 


16 


01357206 PN 


01 03 14.11 


-72 09 14.2 


0.5 


1.34 x 10" 


1 ±6.38 x 10" 


■3 


2049.8 


3.5 x 10" 


13 


17 


01357206 PN 


01 03 37.57 


-72 01 33.2 


0.2 


5.09 x 10" 


1 ±9.62 x 10" 


■3 


17172.3 


2.6 x 10" 


12 


18 


01357206 PN 


01 05 55.38 


-72 03 47.9 


1.6 


3.05 x 10" 


2 ±3.79 x 10" 


-:i 


119.5 


3.9 x 10" 


14 



Notes to column No 1: Observation ID and the instrument with which the detection with the highest likelihood was achieved. 
Notes to columns No 3-5: Position of the XMM-Newton detections with 1 a statistical positional error, except for No 02, for 
which the position of the optical counterpart* is given. 

Notes to column No 6: Count rate for EPIC PN detection (except for entries No 01 and 10 which are EPIC MOS1 detections, 
and t 3 a upper limit for No 02). 

Notes to column No 7: Maximum likelihood (ML) of detection for the total band (0.3 - 10.0 keV). 

Notes to column No 8: Flux (0.3 - 10.0 keV) calculated from the best fit model spectrum, setting the Galactic foreground 
absorption to zero. See Sect . 12.21 for used models. In order to obtain the luminosity, multiply by 4.3 x 10 47 cm 2 . 



(Fig.HSJ and No 17 (Fig.EJ}, the blackbody fit results in 
higher x 2 '- 231.5 for 156 degrees of freedom for source No 

16 and 322.5 for 241 degrees of freedom for source No 

17 (compare to Table |3J). This is because the blackbody 
model fits the low energy tail of the spectrum, but does 
not account for the peaks around 0.6 keV and 0.9 keV, 
which might indicate emission lines from highly ionized 
oxygen and neon, as well as iron lines. For sources which 
were bright enough to obtain a significant spectral fit, the 
results are shown in figures and the model parameters 
yielding the best fit results are listed in the Table to- 
gether with la errors. Moreover, for source No 17 which 



is very bright and shows emission lines most prominently, 
we also used the VMEKAL instead of the MEKAL model. 
In this more elaborate model, the abundance for each of 
the element is a variable fit parameter. We obtained an im- 
proved fit with high oxygen and neon abundances, whereas 
the other elements have values below solar, not differing 
significantly from zero. 

From the spectral models for the emission we were 
able to estimate the flux of the sources in the ROSAT 
band (0.1 - 2.4 keV). The flux was calculated from the 
fitted models, except for four sources which were too faint: 
For No 02, a power law spectrum with T = 0.7 and 
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Table 2. Continued. 



1 


9 


10 


11 


12 


13 


14 


15 


No 


HR1 


HR2 


HR3 


P 

[s] 


do 
["] 


OGLE Name 


[MA93] 


01 


+0.15±0.20 


+0.09±0.18 


-1.00±0.26 


- 


- 


- 


- 


02 








too weak 








03 


+0.41±0.12 


-0.03±0.11 


-0.30±0.17 


59.00±0.02 


1.2 


00545617-7226476 


810 


04 


+1.00±0.23 


-0.15±0.17 


-0.13±0.25 


140.1±0.3* 


_ 


_ 


904 


05 


+0.49±0.05 


-0.29±0.05 


-0.36±0.08 


_ 


1.8 


00571981-7225337 


_ 


06 
06a 


+0.45±0.11 
+0.80±0.31 


+0.04±0.09 
+0.10±0.26 


-0.41±0.11 
-0.54±0.43 




1.9 
3.5 


00573601-7219339 


1020 


07 


+0.19±0.09 


+0.03±0.08 


-0.14±0.10 


281.1±0.2 






1036 


08 


+0.48±0.07 


+0.25±0.05 


-0.06±0.06 


152.34±0.05* 


_ 


_ 


1038 


09 


+1.00±0.33 


+0.44±0.21 


+0.51±0.13 


_ 


4.5 


00581258-7230485 


_ 


10 
10a 


+0.49±0.07 
+0.27±0.09 


-0.15±0.07 
-0.05±0.09 


-0.21±0.10 
+0.12±0.09 


- 


1.2 
1.4 


00592103-7223171 


- 


11 


+0.18±0.14 


-0.19±0.14 


-0.46±0.21 


- 


- 


- 


- 


12 


i i f>n in 1 o 

+1.00±0.13 


n qj-t I i~\ on 

-U.37±0.29 


1 n n 1 1 a co 

+U.01±0.52 








1208 


13 
13a 


+0.37±0.29 
-0.25±0.36 


-0.30±0.25 

n 1 _l n ao 


+0.06±0.41 

n ia I n A 1 








1240 


14 
14a 


+U.14±0.U8 
+0.26±0.12 


1 n 01 i n n*? 

+U.21±0.U7 
+0.06±0.12 


+U.06±0.U7 
-0.10±0.15 


452.2±0.5 
too weak 


2.6 
2.2 


01012064-7211187 


1257 


15 
15a 


+0.17±0.10 
+0.17±0.13 


-0.23±0.10 
-0.36±0.12 


-0.20±0.16 
-0.26±0.30 








12 II 


16 


+0.16±0.07 


+0.11±0.06 


+0.00±0.06 


341.7±0.4 






1367 


17 


+0.28±0.03 


+0.18±0.02 


+0.07±0.02 








1393 


18 


-0.19±0.18 


+0.12±0.18 


+0.22±0.16 








1557 



Notes to columns No 9-11: Hardness ratios as defined in Eq. 
Notes to column No 12: Pulse periods from timing analysis. * New X 
Notes to column No 13: Distance to OGLE object. 
Notes to column No 15: Entry numbers in lMevssonnier fc Azzopardi 



ray binary pulsar! 
1993 



Nu = 0.0 x 10 21 cm" 2 jYokogawa et all l200lh was as- 
sumed to estimate the flux upper limit. For No 01, 09, 
and 12, T = 1.0 and N H = 1.0 x 10 21 cm~ 2 were adopted. 
The resulting luminosity was used to create a long term 
light curve of all the ROSAT and the new XMM-Newton 
data. In the light curves, crosses are used for ROSAT 
PSPC data, triangles for ROSAT HRI data, and dots for 
XMM-Newton EPIC data. Upper limits determined from 
ROSAT observations are plotted as arrows. For the dis- 
tance to the S MC, a mean value of 60 kpc was assumed 
(see review by Ivan den Berghl ^)99) . 



3. Comments on individual HMXBs and 
candidates 

In this section, we present the results on individual 
sources. All sources, which were detected in the four data 
sets and were proven to be HMXBs or candidates, are 
listed in Table |2 

3.1. Source No 1: RX J 005 1.7-7341 

RX J0051.7- 7341 which has been suggested as an XRB 
candidate bv lKahabka et all (|l999j) was only detected in 
MOS1/2 data. In the PN data the source was located on a 
bad column. It is faint, so neither spectral nor timing ana- 
lysis was performed for this so urce. The PSPC count rate 
of 1.64xl0- 3 ±0.68xl0- 3 s" 1 " 



( Kahabka et"ail lT999) dur- 
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Table 2. Continued. 



1 


16 


17 


18 


19 


20 


21 


No 


["] 


HRI 


PSPC 


[HS2000] 


Source ID 


Remarks 


01 


- 


- 


- 


- 


RXJ0051.7-7341 


XRB? 


02 


22.8 




547 


28 


SMCX-2 


HMXB Be, P 


03 


4.5 


058 


241 


31 


RXJ0054.9-7226 


HMXB Be, P 


04 


_ 


_ 


_ 


32 


XMMU J005605. 2-722200 = 2E 0054.4-7237? 


HMXB?, P 


05 


5.6 


_ 


234 


_ 


2E 0055.6-7241, RX J0057.3-7225 


AGN, z = 0.15 


06 
06a 










XMMU J005735. 7-721932 = [Y112000J 19: 


new HMXB r 


07 


8.1 


073 


114 


35 


AXJ0058-720, RXJ0057. 8-7202 


HMXB?, P 


08 


7.3 


074 


136 


36 


RXJ0057.8-7207 


HMXB?, P 


09 


4.3 


076 


_ 


38 


RXJ0058.2-7231 


HMXB Be 


10 
10a 


2.9 
4.7 


081 


218 


- 


RXJ0059.3-7223 


XRB? 


11 


6.6 


088 


123 


- 


RXJ0100. 2-7204 


XRB? AGN? 


12 










XMMU J010030. 2-722035 


new HMac : 


13 
13a 


4.8 
9.8 


093 


132 


42 


RXJ0101.0-7206 


HMXB Be 


14 
14a 


6.3 
6.3 


095 


159 


43 


RXJ0101.3-7211 


HMXB Be, P 


15 
15a 


5.3 
7.4 


096 


Izl 


i i 

44 


Ka J 0101.0— 1 2U4 


TJTi /TVT39 

HMad : 


16 


4.6 


101 


143 


49 


AXJ0103-722, SAX J0103. 2-7209 


HMXB Be, P 


17 


2.6 


105 


106 


50 


RXJ0103.6-7201 


HMXB? 


18 


6.9 




120 


55 


RXJ0105.9-7203 


HMXB? 



Notes to column No 16: Distance to ROSAT source. 

Notes to columns No 17 - 19: Entry numbers in ROSAT HRI catalogue dSasaki et all |200CD . ROSAT PSPC catalogue 

l|Haberl et all l200(t ) . and lHaberl fc Sasakil l|200nl). 

Notes to column No 20: YIT2000: lYokoeawa et"ai] fcOOOl) . 

Notes to column No 21: HMXB: High mass X-ray binary, HMXB?: HMXB candidate, XRB?: X-ray binary candidate, Be: Be 
system, P: pulsar. 



ing the ROSAT observation corresponds to XMM-Newton 
MOS (medium filter) count rate of about 8 x 1CP 3 s _1 . 
This means that the luminosities of the source during the 
ROSAT and XMM-Newton observations were comparable 
(see Tabled- 

3.2. Source No 2: SMCX-2 

SMCX-2 was one of the first t hree X-ray source s 
which were discovered in the SMC llClark et all Il978h . 
It was also de t ected in the HE AO 1 A- 2 experiment 
llMarshall et all Il979l). bu t not in the Einstein IPC 
survey llSeward fc Mitchell Il98lh . In ROSAT obser- 
vations, this transient source was detected only once 



l|Kahabka fc Pietschlll996l) . It is thought to be a Be/XRB, 
since a Be-star was found as its optical counterpart 
i Murdin et aIT.Il979l) . In early 2000, the RXTE All-Sky 
Monitor detected an outburst at the position of SMC X-2 
l|Corbet et all l200ll) and a pulse period of 2 . 374±0.007 s 
was determined l|Corbet fc Marshall l2000l iTorii et all 
l2000l) . 

In the XMM-Newton data (Obs. ID 00842008), there 
was no detection with ML > 10 (see Sect. l2.l|l at the po- 
sition of SMC X-2 which was apparently in low luminos- 
ity state during the XMM-Newton observation. Therefore, 
we performed source detection using the maximum like- 
lihood routine at the position of the optical counterpart 
(SIMBAD): RA = 00 h 54 m 33.4 s , Dec = -73° 41' 04" 



M. Sasaki et al.: XMM-Newton observations of HMXBs in the SMC 
Table 3. Spectral parameters for sufficiently bright sources. 



7 



1 

X 


2 


•i 

o 


4 


K 


a 
u 


7 


u 




q 


1 n 


No 


ID 


r 


N H 
[cm- 2 ] 


kT 

[keV] 


Abund. 
(solar) 


Thermal 
model 


x 2 


dof 


No 


05 


2E 0055.6-7241 


2 59+ - 16 

z ' u3 -0.33 


5.3+o j x 1021 


- 


- 


- 


41.4 


51 


2 


06 


XMMU J005735. 7-721932 


1 42+ - 25 


3.6±i:| x io 21 


- 


- 


- 


84.6 


63 


4 


07 


AXJ0058-720 


1 oi +0 - n 

J - UA -0.U 


3.4±|;| x 10 20 


- 


- 


- 


91.6 


65 


3 


08 


RXJ0057.8-7207 


97+o os 

u - a ' -0.07 


3.0^0 j x 10 21 


- 


- 


- 


161.9 


133 


3 


10 


RXJ0059.3-7223 


1 46+0- 12 


i.8l8:5 x io 21 


- 




- 


117.8 


99 


5 


11 


RXJ0100.2-7204 


2 00+0- 41 


i.7±?:» x io 21 


- 


- 


- 


28.9 


26 


3 


14 


RXJ0101. 3-7211 




3.3l 2 ;t x 10 21 


20+ 09 

u - zu -0.06 


H+ 010 


MEKAL 


82.6 


56 


2 


15 


RX T01 01 6-7904 


i 7 o+0.15 
L - 1 °-0.17 


8 6+ 4 ' 9 x IO 20 
o.u_ 5 6 x lu 








68.0 


58 


5 


16 


AX J0103-722 


i ns+°- 12 
l-U»_o.i 9 


1.9±1;? x 10 21 


n 97+O.O8 
u -^'-0.07 


si+ - 43 


MEKAL 


172.7 


155 


3 


17 


RXJ0103.6-7201 (Model 1) 


n 79 +0.06 
u - ,z -0.07 


1.7t£;° x 10 21 


n 97+O.O3 
u " a '-0.04 




MEKAL 


295.1 


240 


3 




(Model 2) 


n 71 +0.05 


2.3lo'j x 1q21 


n 09+0. 05 
u - oz -0.05 


1.071°;^ (0) 
1.611J:S (Ne) 


VMEKAL 


272.9 


238 


3 



Notes to column No 4: Additional JVh to Galactic foreground iVH.Gai (see Eqs. J5J and @). 
Notes to column No 9: Degrees of freedom. 
Notes to column No 10: Number of used spectra. 



(J2000.0). Since we set the ML limit lower, the source 
was detected with a likelihood of ML = 3.4. The 3 a up- 
per limit count rate obtained from the ML source detec- 
tion routine is 2.33 x IO -3 s _1 . The source counts were 
highest in the B3 band (2.0 - 4.5 keV). In order to es- 
tim ate the flux upper lim it, spectral parameters derived 
bv lYokoeawa et alJ l|200ljl from the ASCA spectrum dur- 
ing the outburst were used: Photon index T = 0.7 for 
a power law spectrum absorbed by a column density of 
iVn < 1.0 x 10 21 cm -2 . This results in an upper limit for 
the un-absorbed flux of 1.5 x 10~ 14 erg cm -2 s _1 , cor- 
responding to L x = 6.5 x 10 33 erg s" 1 (0.3 - 10.0 keV) 
during the XMM-Newton observation in Oct. 2001. 

3.3. Source No 3: RXJ0054. 9-7226 

RXJ0054. 9-7226 is known to be an X-r ay binary pulsar 
with a pulse period of 58.9 69±0.001 s ^Marshall et all 
ll998tlSantangelo et allll998(l and is the only source in our 
samp le, for which the orbi tal period has been measured: 
65 d (Lochner et al ., 1999). In the timing analysis of the 
new XMM-Newton data, the pulse period was verified to 
be 59.00±0.02 s. The folded light curves show variations 
especially above 1.0 keV, and there is no significant change 
in hardness ratios (Fig.[3J) . As can be seen in the long term 
light curve, compared to ROSAT data, the source was ob- 
served in low luminosity state. Due to the low flux, the 
statistics of the spectrum were not high enough and the 
spectrum is thus not discussed here. However, the results 
of the spectral analysis was used to estimate the flux of the 
source (see Table^- The optical counterpart, a Be-star, is 



identified with the v ariable star OGLE 00545617-7226476 
flZebrun et alll200l|) . 



3.4. Source No 4: XMMU J005605. 2-722200 = 
2E 0054.4-7237? 

The error circle of the Einstein source 2E 0054.4-7237 
includes an emission line object. Therefore, it was sug- 
gested as a Be/XRB candidate ([HS2000]). In the XMM- 
Newton data, a source consistent with the position of 
the emission line object was detected (XMMU J005605.2- 
722200) and pulsations from this source was discovered. 
XMMU J005605. 2-722200 is most likely consistent with 
2E 0054.4-7237. The period is 140.1±0.3 s. As can be seen 
in Fig. 01 the pulses in the soft band are narrower than in 
the harder band. 



3.5. Source No 5: 2E 0055.6-7241 

2E 055.6-7241 had been suggested as an XRB candidate 
bv lKahabka et alJ l|l999|) . Timing analysis revealed no pul- 
sations of the X-ray source. Also on longer timescales no 
flux change was verified: The RO SAT PSPC count rat e 
was 7.48 x 10~ 3 ±0.62 x 10~ 3 s" 1 l|Kahabka et allll999|) . 
corresponding to a count rate of 8 x 10 -2 s -1 for XMM- 
Newton EPIC PN (thinl filter). This value is similar to 
the count rate of the XMM-Newton observation, which is 
8.79xl0" 2 ±0.42xl0~ 2 s" 1 . The X-ray spectrum is shown 
in Fig.|SJ It has a photon index of T = 2.591q33, wn i cn is 
higher than for other sources of our sample, and highest 
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00S4:?00i01 F1104 - EPIC PN 





HR1 



Fig. 2. HR1 plotted over HR2 and HR2 over HR3 with 
errors for all the sources in Table |2 except for SMCX-2. 
Circle is used to mark the source No 5 which is an AGN, 
and triangle for source No 11 which is either an XRB or 
an AGN. 



absorbing column density of Ah = 5.31q g x 10 21 cm' 2 
(also see Table |2l . A difference to other sources is also 
seen in the hardness ratios, as the source has a rela- 
tively high HR1 and lower values of HR2 and HR3 (HR1 
= +0.49±0.05, HR2 = -0.29±0.05, HR3 = -0.36±0.08, 
also see Fig.[2J). The high absorption makes HR1 positive, 
whereas HR2 and HR3 are negative due to steeper power 
law spectrum. 

On the DSS2 (red) image, there is a source at the 
X-ray position, which coincid es with the variabl e ob- 
ject OGLE 00571981-72253375 l|Zebrun et alll200l|) with 
B = 19.7 an d R = 17.8 (USNO- A2.0 0150-00625436), i.e. 
B - R : = U. ISasaki et all l)2000|) have shown, that all the 
HMXBs and candidates in the SMC HRI catalogue have 
14 < R < 18 and -2 < B - R < 3, whereas e.g. AGNs 
have R > 16 and B — R> 0. Both the optical magnitudes 
and the X-ray spectra indicate that this source might as 
well be an AGN. Spectroscopy of the optical counterpart 




Pulse Phase 



HRI SMC 058 



50000 

JD 



51000 
2400000 



Fig. 3. Folded light curves and long term light curve of 
RX J0054.9-7226 (source No 3). The hardness ratio is the 
ratio between the count rates in harder band and the count 
rates in softer band fSect. l2.2| . See text for the symbols 
used for the long term light curve. 



bv lDobrzvcki et ahl (|2003) showed that this object is a z 
= 0.15 quasar located behind the SMC. 

3.6. Source No 7: AXJ0058-720 

The pulse period of AX J0058-72 was determined from 
the A SCA data as 280.4±0.3 s l|Yokogawa fc KovamaL 
1998), which we confirmed in the XMM-Newton data: 
281.1±0.2 s. It shows strong pulses in the softer bands and 
its spectrum becomes harder during the 'off' time (Fig.EJ). 
The residuals of the power law fit (Table Eland Fig. El in- 
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0054300101 PN11 - EPIC PN 




II I F1105 - EPIC PN 



Fig. 4. Folded light curves of XMMU J005605. 2-722200 
(source No 4). Hardness ratio as in FigJ3J 



v0>4.::00101 PN20 




0.5 1 2 5 

channel energy (keV) 

Fig. 5. Spectrum of 2E 0055.6-7241 (source No 5). 



dicate the existence of an additional soft component. The 
source has been suggested as a HMXB candidate due to 
the likely optical counterpart, which is an emission line 
object ([HS2000]). 

3.7. Source No 8: RXJ0057. 8-7207 




Pulse Phase 



01 10000201 PN05 




HRI SMC 073 



50000 51000 52000 

JD - 2400000 



RXJ0057.8-7207 is a HMXB candidate with an emis- Fig. 6. Folded light curves, spectrum, and long term light 

sion line object suggested as a likely optical counterpart curve (0- 1 " 2 - 4 keV ) of AXJ0058-720 (source No 7). 

([HS2000]). We discovered pulsations in the new XMM- Hardness ratio and symbols as in FigEl 
Newton data and derived a pulse period of 152.34±0.05 s. 
For this source, a period of 152.098±0.016 s was indepen- 
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cliormel energy (kev) 



HRI SMC 074 



50000 
JD - 2400000 



Fig. 7. Folded light curves, spectrum, and long term light 
curve (0.1 - 2.4 keV) of RX J0057.8-7207 (source No 8). 
Hardness ratio and symbols as in Fig|3 



dently found in Chandra data by iMacomb et al. 
The folded light curves in Fig. show, that there are cor- 
related flux variations in all bands with a significant mi- 
nimum at phase 0.4. Especially in the hard band, there 
is a slow increase and fast decay. Therefore, the hardness 
ratio falls off at phase 0.2 and increases slowly after phase 
0.7. The source spectrum is well reproduced by a power 
law spectrum (see Table |3J) with a significant absorption 
within the SMC or the source itself. As can be seen in the 
long term light curve, there was a weak flare observed by 
ROSAT, whereas the XMM-Newton observation was per- 
formed in a low luminosity state, 5.3 times lower than the 
maximum observed by ROSAT. 

3.8. Source No 9: RXJ0058. 2-7231 

The source corresponding to the optically identified 
HMXB RXJ0058. 2-7231 is very faint, so that no timing 
analysis could be performed. However, the hardness ra- 
tios HRI, HR2, and HR3 indicate, that this source has 
a hard spectrum. Its optical counterpart is a variable Be 
star i n the SMC, OGLE 00581258-7230485 l|Zebrun et all 
l200l|) . From the ROSAT HRI coun t rate of 4.28 x 10" 3 ± 



7Tl l2ooo^ we estimated the cor- 



0.48 xlO-- 3 s" 

responding XMM-Newton EPIC PN (thinl filter) count 
rate: ~ 2 x 10 _1 s _1 . The source was about 3.6 times 
brighter when it was detected by ROSAT than when it 
was observed by XMM-Newton. 



3.9. Source No 10: RXJ0059. 3-7223 

RX J005 9.3-7223 has been su ggested as an XRB candi- 
date bv lKahabka et alJ l|l999h . It was observed by XMM- 
Newton in two pointings. Its spectrum mainly consists of 
a power law component typical for a HMXB with ad- 
ditional features (Fig.[HJ). For this source no pulsations 
were detected. At its posi tion, there is t he va riable star 
OGLE 00592103-7223171 l|Zebrun et all l200l[) . which is 
suggested as the optical counterpart. Its magnitudes are 
B = 17 .4 and R = 14.6 (USNO-A2.0 0150-00660299), 
which gives B — R = 2.8. The R magnitude in particu- 
lar is characteristic for a HMXB (see Sect. l3.5j) . 



3.10. Source No 11: RX JO 100. 2-7204 

At the position of the XMM-Newton detection corre- 
sponding to RXJ0100. 2-7204, a very faint object can be 
found on the DSS2 (red) image. However, there is no 
entry in the USNO-A2.0 catalogue for this source. We 
also looked for information in different catalogues using 
BROWSE of the HEASARC archive, but could not find 
the magnitudes of this optical source. T he X-ray source 
was s uggested as an XRB candidate by iKahabka et alJ 
(1999J). The spectrum of the source is a power law with 
L = 2.00^0 26 an d absorbing column density of Njh = 



1.7 



+0.9 



10 21 cm 2 (Fig. El- Since the probable optical 
counterpart is very faint and the power law photon index 
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Fig. 8. Spectrum and long term light curve (0.1 - 2.4 keV) Fig. 9. Spectrum and long term light curve (0.1 - 2.4 keV) 
of RX J0059.3-7223 (source No 10). Symbols for the long of RX J0100.2-7204 (source No 11). Symbols for the long 



term light curve as in Fig.|3| 



term light curve as in Fig.|3J 



is higher than for most of the other sources presented here, 
it can not be ruled out that this source is an AGN (also 
see Table EJ. 

3.11. Source No 13: RXJ0101. 0-7206 

The Be/X-ray binary RXJ0101. 0-7206 showed a luminos- 
ity of - 3x 10 33 erg s" 1 in the ROSAT band (0.1 - 2.4 keV) 
during two XMM-Newton observations. It was about 60 
times fainter than at the maximum observed by ROSAT 
(Fig.llOJl. Pulsations with a p eriod of 304.49±0.13 s were 
discovered in Chandra data ijMacomb et all [2003^ . This 
period could not be verified in the X MM-Newton observa- 
tion, because the source was too faint. lEdee fc Codl|2003|) 
presented results on the optical analysis of likely counter- 
parts, discussing two objects (No 1 and 4) in the ROSAT 
PSPC error circle. They conclude that the optical coun- 
terpart is object No 1 which is confirmed to be a Be star. 
This object is also the only optical source, which can be 
found on the DSS image within the XMM-Newton 1 a 
error circle. 



3.12. Source No 14: RX JOIOI. 3-7211 

The ROSAT source RX J0101. 3-7211 is the first X-ray 
binary pulsar of which the discovery was based on XMM- 



HRI SMC 093 



50000 
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2400000 



Fig. 10. Long term light curve (0.1 - 2.4 keV) of 
RXJ0101. 0-7206 (source No 13). Symbols as in FigJSI 

Newton data. It was covered in two additional observa- 
tions finding the source again in a low inten sity state. The 
pulse period of 455±2 s llSasaki et all boOD was verified in 
the new data of the observation ID 01100002: 452.2±0.5 s. 
During the observation ID 01357206, the source was too 
faint for a timing analysis. The folded light curves show 
strong variation in all bands (Fig. lllll . The spectrum of 
the source becomes harder during pulse minimum. The 
spectrum is well fitted with a soft thermal component de- 
scribed by a MEKAL model (kT = 0.20±^ keV) with a 
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Fig. 11. Folded light curves, spectra, and long term light 
curve (0.1 - 2.4 keV) of RX J0101. 3-7211 (source No 14). 
Hardness ratio and symbols as in FigEl For the spectra, 
solid lines are used for the data of the obs. ID 01100002, 
and dashed lines for obs. ID 01357206. 
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Fig. 12. Spectra and long term light curve (0.1 - 2.4 keV) 
of RXJ0101. 6-7204 (source No 15). For the spectra, solid 
lines are used for the data of the obs. ID 01357206, and 
dashed lines for obs. ID 01100002. Symbols for the long 
term light curve as in Fig. El 



low metal abundance (0.11^0 u times solar) and a power 
law component absorbed by a high column density (Table 
EJ). The optical counterpart (OGLE 01012064-7211187) is 
a Be-star. 

3.13. Source No 15: RX J0101. 6-7204 

The Be/XRB candidate RXJ0101. 6-7204 with an emis- 
sion line star at the ROSAT PSPC and HRI positions 
([HS2000]), was observed in two XMM-Newton pointings. 
Its spectrum and long term light curve are shown in 
Fig-El The spectrum can be modelled as a moderately 
absorbed power law. No pulsations were discovered. 

3.14. Source No 16: AXJ0103-722 

For the Be/X-ray binary AX J O 103-722 a pulse p eriod of 
345.2±0.1 s was determined bv llsrael et alJ l|l998|) . In the 
XMM-Newton data, pulsations were confirmed with a pe- 
riod of 341.7±0.4 s. The folded light curves show strong 
variation below 2.0 keV fFig. H3|) . whereas in the hard 
band, the variations are strongly reduced. The spectrum 
is well reproduced with a power law and a thermal com- 
ponent (see Table OJl. The MEKAL model for the thermal 
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Fig. 13. Folded light curves, spectrum, and long term 
light curve (0.1 - 2.4 keV) of AXJ0103-722 (source No 
16). Hardness ratio and symbols as in FigGH 
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Fig. 14. Spectrum and long term light curve (0.1 - 
2.4 keV) of RX J0103.6-7201 (source No 17). Symbols for 
the long term light curve as in Fig.|31 



component yields kT = 0.27 



+0.08 
0.07 



keV and metal abun- 



dances of 0.31^0.17 with respect to solar. 



3.15. Source No 17: RXJ0103.6-7201 

For the HMXB candidate RX J0103. 6-7201 ([HS2000]), 
an acceptable fit was obtained for the spectrum with a 
power law and a thermal component (Table EJ. Modelling 
the thermal component with MEKAL, we obtained kT = 
0.271q'q4 and metal abundances of 0.77lg'2g times solar. 
Since the source was bright, we also used the VMEKAL 
model instead of the MEKAL model, which allows to de- 
termine the abundance for each of the elements. This re- 
sulted in an improvement of the fit, the model reproducing 
the peaks around 0.6 and 0.9 keV. The photon indices T 
and the absorbing column densities Ah for both fits are 
comparable, as can be seen in Table El Also the temper- 
ature values kT agree well for MEKAL and VMEKAL 
within the 1 a errors. The spectrum with the power law 
+ VMEKAL fit is shown in Fig.^1 The comparison to 
ROSAT data shows that this source was in high lumi- 
nosity state during the XMM-Newton observation with 
L x = 1.1 x 10 36 erg s" 1 (0.3 - 10.0 keV). In spite of the 
high photon statistics with 3,300 counts, no pulsations 
were discovered. Also the analysis of the events separated 
into soft, medium, and hard band revealed no pulsations. 
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Fig. 15. Spectrum of XMMU J005735. 7-721932 (source 
No 6). 



3.16. Source No 18: RXJ0105. 9-7203 

RXJ0105. 9-7203 is a HMXB candidate, coinciding with 
an emission line object. Since the source was very faint 
during the XMM-Newton observation, the photon statis- 
tics are very low and no timing analysis was possible. The 
PSPC count rate derived from the ROSAT observatio n 
was 4.01 x 10" 3 ± 0.56 x 10" 3 s" 1 l|Haberl et all l2000j) . 
corresponding to - 5 x 10 -2 s _1 for XMM-Newton EPIC 
PN (thinl filter). With a count rate of 3.05 x 10~ 2 ± 
0.28 x 10~ 2 s _1 (Table |2J), the source was fainter during 
the XMM-Newton observation. 



4. Sources No 6 and 12: New HMXB candidates 

To identify a HMXB, it is crucial to find an optical coun- 
terpart and confirm that it is an early-type star. If an 
emission line object is found at the position of a hard X- 
ray source and other objects are ruled out as counterpart, 
the source is presumably a Be/XRB. Cross-correlating the 
XMM-N ewton source list with the emissio n line star cat- 
alogue of lMevssonnier k. Azzopardl l)l993|) . we discovered 
two new sources which met these criteria. 

XMMUJ005735. 7-721932 (source No 6) is found 
at the position of [MA93110 20 and likely c o incide s 
with the source No 19 in lYokogawa et all l|2000|) . 
XMMU J010030.2-722035 (source No 12) is associated 
with [MA93]1208. Both X-ray sources are very faint. Only 
for XMMUJ005735.7-721932, we had enough counts to 
extract a spectrum (Fig.ll5|l. The best fit model is a 
moderately absorbed power law (Table |3J) . Furthermore, 
Chandra data showed that this so urce has pulsed emis - 
sion with a period of 564.81±0.41 s llMacomb et alll2003|) . 
This period could not be confirmed in the XMM-Newton 
data. 

As indicated by the hardness ratio HR1, these sources 
have a hard spectrum. Therefore, these two sources are 
suggested as new Be/XRB candidates. For further inves- 
tigation, we need to perform follow-up observation in the 
optical band in order to verify if the emission line objects 
are Be stars. 
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Fig. 16. Histogram of luminosities [erg s x ] of the 
HMXBs and candidates in the XMM-Newton band (0.3 
- 10.0 keV). The upper limit for SMCX-2 is shown with 
dashed line. 

5. Discussion 

The comparison between the XMM-Newton sources de- 
tected with ML > 10 and other X-ray catalogues 
llWarig fc Wil 119921 iKahabka et all Il999t lHaberl et all 
1200(1: ISasaki et all l200ff demonstrates that we detected 
all known HMXBs and candidates which exist in the four 
observed fields, except for SMCX-2 which was very faint. 
SMC X-2 was marginally detected with a likelihood of ML 
= 3.4, and we derived an upper limit of 6.5 x 10 33 erg s _1 
(0.3 - 10.0 keV). The luminosities of all the other sources 
in the 0.3 - 10.0 keV band are high er than 8 x 10 33 erg s" 1 
at an assumed distance of 60 kpc Ijvan den Berg 
as is shown in Fig. 1161 As we have seen in the long term 
light curves, most of the sources were in quiescence during 
the XMM-Newton observations, whereas they were mostly 
detected during outburst by previous missions. This indi- 
cates that all known HMXBs in the SMC have luminosities 
higher than ~ 7 x 10 33 erg s^ 1 in quiescence and can be 
detected by XMM-Newton in observations with an expo- 
sure of about 15 ks. Consequently, we have an extensive 
set of HMXBs for studying their properties. 

In order to visualise the spectral characteristics of the 
HMXBs, we plotted the hardness ratios HR1, HR2, and 
HR3 in Fig. |2 The high absorption in XRBs causes pos- 
itive values for HR1, while HR2 and HR3 have small ab- 
solute values around zero. AGNs typically show steeper 
X-ray spectra than HMXBs. Therefore, the two source 
classes can be disentangled using hardness ratios. This 
can also be applied to classification work on other nearby 
galaxies. 

5.1. Pulsations and soft energy emission 

X-ray spectra of high mass X-ray binaries below 10 keV 
can be in general modelled as a power law with a pho- 
ton index of T = — 2. For HMXBs located far away 
from the Galactic plane or in the Magellanic Clouds, 
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the interstellar absorption in the line of sight is low, 
and an additional soft spectral component was discov- 
ere d in supergian t X-ray binary sy s tems like SMC X- 
1 {Marshall et all Il983t IWoo et all Il995h or LMCX 



4 I Woo et all Il996ft . as well as in Be/X-ray bin ary 
systems like RXJ0 059. 2-7138 llKohno et all 12000ft or 
EXO 053109-6609.2 llHaberl et alll2003h . 

In our sample of HMXBs and candidates in the SMC, 
pulsations were confirmed for six sources. Studying the 
pulsations and hardness ratio changes in different bands, 
we found that there are different types of pulsations. 
Furthermore, four out of these six were bright enough to 
allow us to test the existence of a soft component in their 
spectra. The pulsating sources of our sample can be di- 
vided into four groups: 

1. There are pulsations in all bands and the ratios be- 
tween the harder and softer bands are almost con- 
stant (RX J0054.9-7226, Fig.E|and XMMU J005605.2- 
722200, Fig.EJ. Because of low photon statistics, spec- 
tral analysis of these sources yielded no significant re- 
sults. 

2. Pulsations are discovered below 2 keV (AX J0058-720, 
Fig.Eland AX J0103-722, Fig.Q]|. There might be pul- 
sations also in the hardest band, although not signif- 
icant due to low statistics. The hardness ratios seem 
to become higher during pulse minimum. The spec- 
tra of these sources include a thermal component. If it 
is confirmed that pulsations are in fact existent only 
in the soft band, this will suggest that their origin is 
of small size, probably locally illuminated surface or 
surroundings of the neutron star. 

3. In the case of RXJ0057. 8-7207 (Fig.0), there are flux 
variations in all bands. Most pronounced pulsations 
are found above 2 keV with correlated increase of the 
hardness ratio, which follows the pulse shape of the 
hard band. Since the spectrum is a perfect power law, 
the coincidence of the maximum of the hardness ratio 
and that of the pulse might indicate a variation of the 
absorption (Ah). 

4. RXJ0101. 3-7211 (Fig-inj shows clear pulsations in all 
bands and becomes harder significantly at pulse mini- 
mum. In its spectrum, there is a low energy excess 
besides the power law component, which can be mo- 
delled as a thermal emission. The pulsations in the soft 
band might be caused by Ah variations as well as by 
changes in the soft emission component. 

The low energy component, which seems to be thermal 
was also found in the spectrum of RXJ0103. 6-7201. This 
source is variable on long timescales, as can be seen in 
Fig-EJ However, pulsations on short timescales were not 
discovered, although the source was brighter during the 
XMM-Newton observation than in former observations. 



5.2. Origin of the soft emission 

The low energy component in the spectrum of th e su- 
pergiant systems SMCX-1 {Marshall et all Il983ft and 



LMC X-4 {Woo et all Il996ft was modelled as blackbody 
emission or thermal Bremsstrahlung which arises from the 
stellar wind of the supergiant, the accretion disk, or the 
fan-beam of the accr etion column clos e to the neutron 
star surface. However. IPaul et all {2002) pointed out that 
a power law nature is most probable for the soft emission. 
They derived that the pulse shape of the soft emission 
from SMC X-l is sinuso idal, similar to the soft ene rgy light 
curve of Her X-l (e.g. lOosterbroek et ail Il997h . In our 
Galaxy, the supergiant system Vela X-l is thought to show 
emission from the atmosphere and stellar wind of the com- 
pani on as well as from the ga s stream towards the neutron 
star {Haberl fc Whitel . ll990t and references therein). High 
resolution spectroscop y of Galactic HMXBs lik e CenX-3 
(with Chandra HETG, IWoidowski et all 12002ft or Her X- 
1 (with XMM-Newton RGS. IJimenez-Garate et all l2002j) 
resolved fluorescent lines and hydrogen- and helium-like 
lines of elements from Ne to Fe. The line fluxes of Cen X-3 
are consistent with recombination radiation from photo- 
ionised and collisionally ionised plasma as well as resonant 
line s cattering in photo-ionised plasma I Woidowski et alj . 

As for the Be/X-ray binary systems, iKohno et alJ 
[|2000ft analysed both ASCA and ROSAT data of 
RX J0059. 2-7138 and found that there is a soft compo- 
nent in the spectrum, which can be modelled as a thermal 
emission with kT = 0.37 keV. Below 2.0 keV, the source 
shows no pulsations. Therefore they argue that the soft 
emission originates from a large region comparable to the 
full binary system. Using an XM M-Newton observa tion 
of a northern field in the LMC, lHaberl et all {2003ft ex- 
tracted emission from the Be/X-ray binary EXO 053109- 
6609.2 and showed that there are strong pulsations above 
0.4 keV. In the spectrum there is a low energy thermal 
component, which is believed to arise from the equatorial 
disk around the Be star, illuminated by the X-ray source. 

The origin of the soft emission from HMXBs is not 
clearly understood. One would expect that there are dif- 
ferences between a supergiant and a Be system. Most of 
the HMXBs which have been studied in detail (since they 
are located in the Milky Way and therefore closer) are 
supergiant systems, whereas the sources in the SMC we 
are confronted with, are Be systems. In Bc/XRBs, the 
neutron star and the Be star are thought to form a bi- 
nary system with an extended orbit. This makes the stel- 
lar material in the equatorial disk around the Be star as 
the origin of the soft pulsed emission rather implausible. 
The HMXBs in the MCs are ideal objects to study the soft 
part of their spectrum, since the absorption by Galactic 
foreground matter is low in the direction of the MCs. The 
existence of a soft thermal component in the spectrum and 
pulsations below 1-2 keV in our data indicates that the 
size of the origin of the soft emission is not as large as is as- 
sumed for e.g. RX J0059. 2-7138. In addition to timescales 
and luminosities, a crucial parameter for the physical pro- 
cesses responsible for this emission is the magnetic field 
of the neutron star. In order to clarify the conditions in 
which the soft component is produced, at least we need 
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to get information about the orbital motion and about a 
possible orbital phase dependence of the total source spec- 
trum as well as the pulsed emission. As for the SMC Be 
systems discussed here, the orbital period is known only 
for one source. 

5.3. OB systems vs. Be systems 

In the last few years, the number of known Be/XRBs in 
the SMC increased drastically based on temporal studies 
of hard X-ray sources and optical observations. In order to 
identify an X-ray source as a HMXB and clarify the nature 
of the mass donor star, we need to perform spectroscopy 
of the optical counterpart. Since most of the HMXB can- 
didates which are known now are correlated to emission 
line objects, we expect that additional Be/XRBs will be 
found in the near future. This will further increase the ra- 
tio between the Be systems and the OB systems among 
the HMXBs in the SMC. Be/XRBs are thought to evolve 
from binary systems in about 1.5 x 10 7 yrs, whereas su- 
pergiant systems evolve faster due to the high mass of 
the companion star. The large number of Be/XRBs sets 
constraints on the secondary star formation in the SMC, 
making a burst some 10 7 yrs ago most likely. 

6. Summary 

We analysed XMM-Newton EPIC PN and MOS 1 /2 data 
of four pointings towards the SMC. One observation co- 
vered the field around the HMXB SMC X-2 in the south, 
whereas the fields of view of the other three are located 
in the northern part of the main body of the SMC. In 
total, there were 15 detections which were identified as 
known HMXBs or XRB candidates. For SMC X-2 which 
was faint during the observation, a flux upper limit of 
1.5 x 1CT 14 erg cm" 2 s" 1 (0.3 - 10.0 keV) was derived. 
We found two new sources (XMMU J005735. 7-721932 
and XMMU J010030.2-722035) which have a hard spec- 
trum and positionally coincide w ith emission line objects 
l|Mevssonnier fe Azzopardi [l99^l . These sources are pro- 
posed as new HMXB candidates, probably Be systems. 

Four sources in our list were known to show pulsed 
emission and pulse periods had been determined in for- 
mer observations. In this work, the pulse periods were 
confirmed for all four sources. Furthermore, we discovered 
that two other sources which had been proposed to be 
Be/XRB candidates, show pulsations: XMMU J005605.2- 
722200 with a pulse period of 140.1±0.3 s and 
RX J0057.8-7207 with 152.34±0.05 s. 

Spectral analysis of the sources was performed. For 
faint sources, a good fit was obtained with a single power 
law spectrum. However, for three brighter sources, we 
could show that there is a significant low energy excess in 
the spectrum, if we only assume a power law. The spectra 
indicate emission line features, suggesting that the emis- 
sion is thermal. This soft component was modelled as ther- 
mal emission, yielding temperatures of 0.2 - 0.3 keV. The 
abundances in the emitting plasma are below solar values, 



but c omparable to typical SMC values l|Russell fe DopitaL 
Il992ft : for RX J0101. 3-7211 it is 0.11+^? times solar, and 
for AX J0103-722 best fit is obtained with 0.31+0.17 times 
solar. The errors are 1 a values. Only for RXJ0103.6- 
7201 the abundance is higher with 0.77+o' 2 g with respect 
to solar. 

The flux of the sources in the MCs is low compared 
to the bright (L x = 10 37 " 38 erg s _1 ) HMXBs in our 
Galaxy, making it difficult to perform a detailed analy- 
sis of their soft emission. However, the sources in the MCs 
have the advantage of low Galactic absorption. This al- 
lows us to study the thermal emission from a large sample 
of HMXBs and to increase the understanding of the in- 
teraction between X-rays from the compact object and 
the ambient stellar matter. It is also important to verify 
if there is a change in temperature or emissivity, which 
is related to the orbital phase of the binary system. Due 
to the improved time resolution and sensitivity, there is a 
large detection potential for new pulsating XRBs in fur- 
ther XMM-Newton observations. 
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